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A relation between a class of stationary points of the energy landscape of continuous spin mod- 
els on a lattice and the configurations of a Ising model defined on the same lattice suggests an 
approximate expression for the microcanonical density of states. Based on this approximation we 
| conjecture that if a 0(n) model with ferromagnetic interactions on a lattice has a phase transition, 

its critical energy density is equal to that of the n = 1 case, i.e., a system of Ising spins with the 
same interactions. The conjecture holds true in the case of long-range interactions. For nearest- 
neighbor interactions, numerical results are consistent with the conjecture for n = 2 and n = 3 in 
<r | t three dimensions. For n = 2 in two dimensions (XY model) the conjecture yields a prediction for 

the critical energy of the Berezinskij-Kosterlitz-Thouless transition, which would be equal to that 
of the two-dimensional Ising model. We discuss available numerical data in this respect. 
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The stationary points of a function of many variables are the points of vanishing gradient and play a relevant role 
in quite a few theoretical methods in physics. When the function is the energy of a many-body system these methods 
are referred to as "energy landscape methods" [ll. Examples of applications include clusters [l| , disordered systems 
and glasses Q , biomolecules and protein folding . Energy landscape methods allow to estimate dynamic as well as 
static properties. As far as equilibrium is concerned, the classic application is Stillinger and Weber's thermodynamic 
formalism QJ, where minima are the stationary points to be considered. Later, all the stationary points of the energy, 
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including saddles of any index 31 [ have been taken into account, for instance to characterize glassy behavior It 
was further realized that stationary points of the Hamiltonian are connected with topology changes of the phase 
space accessible to the system, leading to the conjecture that some of them are at the origin of thermodynamic 
phase transitions d, Q ; quite some research activity followed, reviewed in d, Q . Although in equilibrium statistical 
mechanics phase transitions in systems with non-fluctuating particle numbers have been mainly studied within the 
canonical ensemble, the relation between stationary points of the Hamiltonian and equilibrium statistical properties 
is more transparent in a microcanonical setting. This can be understood in an intuitive way by observing that the 
entropy is defined as i 32j s(e) — logw(e), where e = E/N is the energy density and uj is the density of states. For a 
system with N degrees of freedom described by continuous variables the latter can be written as 

where T is phase space and dT its volume measure, E £ is the hypersurface of constant energy E — Ns, and c?S 
stands for the A" — 1-dimensional Hausdorff measure. The rightmost integral stems from a coarea formula [l0(. At a 
stationary point, VH = and the integrand diverges, so that its contribution to U) is clearly important. Indeed, the 
density of states is nonanalytic at stationary values of the energy for any finite N, and so is the entropy, at variance 
with the canonical free energy which may develop nonanalyticities only in the thermodynamic limit A" — > oo [Tl| . 
Microcanonical nonanalyticities at finite N are in one-to-one correspondence with stationary configurations; however, 
the "strength" of such nonanalyticities generically decreases linearly with N, i.e., the first k derivatives of the entropy 
are continuous, where k is O(N) (KSS theorem [12, [Uj])- The usual thermodynamic functions are given by low-order 
derivatives of the entropy, so that these nonanalyticities can be observed only for very small N from noisy data. In 
the thermodynamic limit most of these nonanalyticities disappear. Only those singularities (if any) that survive (or 
appear) in that limit are typically associated to thermodynamic phase transitions and coincide with the canonical 
nonanalyticities if equivalence of statistical ensembles holds. This may suggest that finite-A^ nonanalyticities are 
totally unrelated to thermodynamic phase transitions. However, this is not true, since Franzosi and Pettini showed 
that stationary points are a necessary condition for phase transitions to occur, at least for systems with short-range 



interactions [l4 |. Moreover, a possible scenario allowing some finite- A" singulartities to survive in the thermodynamic 
limit has been depicted in [12J. Therefore, the relation between the energy landscape, i.e., the stationary points of 
the Hamiltonian, and thermodynamic phase transitions is still an open problem. 

The purpose of the present paper is to apply an energy landscape analysis to 0(n) spin models. Loosely speaking, 
we take seriously the idea that the stationary points contribution to the density of states is the most important one, 
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and by considering a particular class of stationary points we construct an approximate form of the density of states 
which suggests a relation between the phase transitions occurring in these models. Let us consider a classical isotropic 
spin model defined on a lattice (or more generally on a graph) with Hamiltonian 

N N n 

« (n) = - E J a s i ■ s i = - E J ^ E s i s " > ( 2 ) 

— 1 2)J=1 Q — 1 

where i and j run over the iV lattice sites and the classical spin vectors Si — (S^ 1 , . . . , S") have unitary norm, 
i.e., 53o=i (S?) 2 = 1 Vi = 1, ...,N. The real matrix Jy dictates the interactions; in case they are long-ranged 
a normalization is understood such as to obtain an extensive energy, using e.g. the Kac prescription |15| . The 
Hamiltonian ^ is globally invariant under the 0(n) group. In the special cases n = 1, n = 2, and n = 3, one obtains 
the Ising, XY, and Heisenberg models, respectively. The case n = 1 is even more special because O(l) = Z2 is a 
discrete symmetry group. In this special case the Hamiltonian © becomes the Ising Hamiltonian 

N 

n {1} = - E Jaw . ( 3 ) 

where <Xj = ±1 Vi. In all the other cases n > 2 the 0(n) group is continuous; each spin vector Si lives on an n — 1 unit 
sphere §" _1 . Let us now consider the stationary configurations of for n > 2, i.e., the solutions S = (Si, . . . , 5jv) 
of the N vector equations VH^V = 0. The latter can be written as nN scalar equations, 

N 

-E JkjSf + XkSZ = , o = l,...,n, fc = l,...,JV, (4) 

3=1 

where the A's are iV Lagrange multipliers, plus the N nonlinear constraints (Sf) 2 = 1, which prevent the 

above equations from being easily solved. However, a particular class of solutions can be found by assuming that 
all the spins are parallel or antiparallel: Sj = ■ ■ ■ = S,™ -1 = Vi. In this case, the N(n — 1) equations (0| with 
a = 1, . . . , n — 1, corresponding to the first n — 1 components of the spins, are trivially satisfied. As to the n-th 
component, the constraints (S 1 ") 2 = 1 imply S™ — <Ti Vi, so that the remaining N equations read as 

N 

~E J *i a i + Afe(Tfe =0 ' * = 1,..-,JV- (5) 

4=1 

The above equations are satisfied by any of the 2 N possible choices of the er's provided one puts Afe = 
(^2^ =1 JkjCj) /vk, k = 1,...,N. The Hamiltonian ^ becomes the Ising Hamiltonian ^ when the spins belong 
to this class of stationary configurations. Therefore we have a one-to-one correspondence between a class of station- 
ary configurations of the Hamiltonian ([2} of a 0(n) spin model and all the configurations of the Ising model ([3]), 
i.e., the Ising model defined on the same graph with the same interaction matrix J^; the corresponding stationary 
values are just the energy levels of this Ising Hamiltonian. We shall refer to the class of stationary configurations 
Si = (0, . . . , 0, <7j) Vi = 1,...,N as "Ising stationary configurations". There will be also other stationary configu- 
rations; nonetheless, the 2 N Ising ones are a non-negligible fraction of the whole, especially at large ./V because the 
number of stationary points of a generic function of N variables is expected to be 0(e ) [la ]. 

The above results hold for 0(n) and Ising models defined on any graph. From now on we shall restrict to regular 
rf-dimensional hypercubic lattices and to ferromagnetic interactions > 0. In this case, in the thermodynamic limit 
N — » 00 the energy density levels of the Ising Hamiltonian (J3|), H^ 1 ' (01, a^)/N Vcr^ = ±1, become dense and 
cover the whole energy density range of all the 0(n) models. This suggests that Ising stationary configurations are 
the most important ones, so that we may approximate the density of states o/™)(e) of an 0(n) model in terms of 
these configurations. To this end, let us first rewrite Eq. (fTJ as 
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where p runs over the 2 N Ising stationary configurations and U p is a neighborhood of the p-th Ising configuration 
such that {Up} 2 =1 is a proper partition of the configuration space F = (S" -1 )^, that coincides with phase space for 
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spin models ([2]). Since Ising configurations are isolated points in the configuration space of a 0(n) model, such a 
partition always exists. Let us now introduce two assumptions allowing to write Eq. (J6)) in a more transparent, albeit 
approximate, way. (i) At a given value of e, the largest contribution to wW(e) is likely to come from those U p such 
that nt n) (p) = Ne, because if ~UW{q) ^ Ne then \VU^ n \x)\ ^ \/x e U q n E e , unless a zero in \VH^(x)\ comes 
from a stationary configuration which does not belong to the Ising class. According to our previous considerations, 
we assume that non-Ising stationary configurations can be neglected. We shall therefore consider only stationary 
configurations at energy density e in the sum (|BJ). (ii) We shall assume that the integrals in Eq. ([6]) depend only on 
e, i.e., the neighborhoods U can be deformed such as 

1 ' K = g^(e) (7) 



'u P n^ |VW(»>| J Uq nv c \VU(")\ 
for any p, q such that 1-L <Jl \p) = 'H^ n \q) = Ne. Hence, using assumptions (?) and (ii), Eq. becomes 

uW (e) ~ g (n) (e) ^ <5 (p) - We] ■ (8) 
p 

The sum on the r.h.s. of Eq. (J5J is over Ising configurations, so that it equals the density of states of the corresponding 
Ising model, that we shall denote by u/ 1 -*^). We can thus write 

Lu^(e)^uj (1) {e)g in) (e). (9) 

Were Eq. © exact, it would imply that if tu^^e) is nonanalytic at e = e c , then also uj^ n \e) is nonanalytic at e = e c 
for any n, unless the function g^ n \e) precisely cancels this nonanalyticity, which seems a rather special case. We do 
not expect Eq. ((9]) to be exact, even in the thermodynamic limit N — > oo, unless, again, g{e) has some very special 
features: with a generic g{e) a density of states of the form (|9]) would not reproduce the known critical exponents 
of the 0(n) universality classes [l7|. However, it can be shown that with a generic g(e) Eq. (j9|) correctly implies a 
negative value for the specific heat critical exponent of O(n) spin models (i.e., the specific heat of continuous models 



does not diverge at criticality, but rather has a cusp-like behavior). This is a common feature of 0(n) models [17 1 
and reinforces the belief that the approximation ([9]) , although rather crude, may properly capture the main features 
of the nonanalyticities of the density of states when N — > oo, as the location of such nonanalyticies. Therefore we put 
forward the following 

Conjecture. If a 0(n) spin model defined on a d- dimensional hypercubic lattice with Hamiltonian @) and ferromag- 
netic interaction matrix Jij > has a phase transition, its critical energy density e c — E c /N is equal to that of the 
n = 1 case, i.e., a system of Ising spins with the same interactions. 

The above conjecture concerns the critical value of the control parameter of the microcanonical ensemble, the energy 
density, and says nothing about critical temperatures, which may well be different — and typically are — at different 
n. 

We now discuss known results, both analytical and numerical, in order to assess the validity of this conjecture in 
some particular cases. The results we were able to collect are reported in Table |U The conjecture is true for for 
systems with long-range interactions on d-dimensional lattices, Jij — N^ a ' d ' 1 \i — j\~ a with 0<a<d;a = 0is 
the mean-field case of models defined on complete graphs with the same interaction strength between any two sites, 
Jij = 1/N. All these systems have a mean-field-like phase transition at the maximum value of e (s c = with our 



choice of units), with critical temperatures T c = 1/n [181 ]. We stress again that critical energy densities are equal 



but critical temperatures are not and depend on n [24|. As to systems with nearest- neighbor interactions, the energy 
density range is e G [—d, d] with our choice of units. The conjecture is true for d = 1 at any n, although this case is 
somehow trivial because there is no transition at finite temperature. For d = 2, the Mermin- Wagner theorem rules out 
a long-range-ordered phase for any n > 1. However, a remarkable transition between a disordered and a quasi-ordered 
phase occurs for n — 2 (XY model) , usually referred to as the Berezinskij-Kosterlitz-Thouless (BKT) transition [25[ . 
In Table J] we report the best recent estimate of the critical temperature obtained by Hasenbusch and coworkers (see 
e.g. [2(| and references quoted therein) and the corresponding critical energy density (estimated from a MonteCarlo 
simulation of a system with 256 x 256 spins [3]). The difference between this value and the exact value of the critical 
energy density of the Ising model on a square lattice is around 2%. This difference, though small, appears significant 
since it is orders of magnitude larger than the statistical error on the numerical estimate of the energy. Based on this 
result one should conclude that the conjecture is not verified in the case of the XY model in d = 2. However, we 
are comparing an exact result in the thermodynamic limit with a numerical estimate of the energy on a finite lattice, 
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TABLE I: Comparison of critical energy densities e c and critical temperatures T c for ferromagnetic models with long-range 
(LR) interactions (first row) and nearest-neighbor interactions on a d-dimensional hypercubic lattice (all the other rows). 
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2.269. . . 
0.8929(1) 


exact solution 6 
numerical [19, 20] 






Ising 


-0.991(1) 


4.5112(3) 


numerical [21] 


d = 3 




0(2) 


-0.991884(6) 


2.2016(7) 


numerical [22] 






0(3) 


-0.9896(1) 


1.44298(2) 


numerical [23] 



"Long-range interactions Jij tx \i — j\ a with < a < d. 
'Exact values are e c = — v/2 and T c = — ^ — 



whose statistical accuracy does not consider the systematic error due to the finite size effects, which could be quite 
large in this particular case [2(| [27]]. Moreover, also the precise determination of the critical temperature of the BKT 
transition is a subtle and difficult task due to its elusive nature. This is witnessed by the remarkable spread of values 
of T c reported in different papers: the summary given in Ref. |27j shows that estimated critical temperatures vary 
in the interval [0.88,0.99] while Ref. [28| gave [0.85,0.95] as confidence interval for T c . The energy values given in 
Ref. [l9| corresponding to both these temperature intervals do contain the Ising value e c = — V2; for instance, the 
temperature interval [0.85,0.95] corresponds to e c <G [—1.48,-1.38]. We thus believe that the available data are not 
conclusive as far as a confirmation of the conjecture is concerned in this particular case. 

For d — 3 the comparison is entirely between simulation outcomes, since no exact solution exists even for the Ising 
case. Results in Table Q] show that the critical energy measured for a 0(2) spin system (XY model) [22| is clearly 
consistent with that measured for the Ising case [21|. The difference between the estimated e c of the 0(3) case 
(Heisenberg model) [23[ and that of the Ising model is less than 1.5 times the error on the latter. Therefore the two 
estimates are consistent if one considers quoted errors as standard statistical errors. 

In conclusion, by considering a special class of stationary points of the energy landscape of 0(n) spin models we 
have proposed an approximate form of the microcanonical density of states and conjectured that the critical energy 
densities of 0(n) models with n > 2 equal those of the corresponding Ising models. Available analytical and numerical 
data are consistent with the conjecture, with the exception of the 0(2) case in d = 2. The latter, however, seems the 
most interesting one. On the one hand, we already noted that the best recent estimates suggest the conjecture does 
not hold in this case but available data seem not conclusive. On the other hand, in this case the conjecture yields an 
exact prediction for the critical energy of the BKT transition in a model that is not exactly solvable. Therefore, in 
our opinion, more precise numerical estimates of the critical energy of the BKT transition would be very interesting, 
as well as estimates of other critical parameters made using the conjectured value of e c . Were the validity of the 
conjecture ruled out in the XY case in d = 2, it would confirm once more the special nature of the BKT transition 
and our conjecture might still be valid in a weaker form, i.e., restricted to 0(n) spin models with symmetry-breaking 
ferromagnetic phase transitions. Conversely, besides yielding an exact value for the critical energy of the BKT 
transition, a confirmation of the conjecture for the XY model in two dimensions might hint at a stronger version of 
the conjecture. Such a stronger conjecture would be that any 0(n) spin model on a d-dimensional hypercubic lattice 
has a phase transition precisely at the critical energy density of the d-dimensional Ising model, even for n > 2 in 
d = 2. We note that the presence of a phase transition in 0(n) models at any n in two dimensions has already been 
suggested by Patrascioiu et al. [29}. Although, to the best of our knowledge, no direct evidence of such a transition 
has been found yet, the possibility remains that the transition exists but is weak and elusive: our conjecture might 
help in finding it, suggesting where to look at. Finally we note that our conjecture might be true for a completely 
different reason than the one argued in the present paper: the phase transitions in 0(n) models may all come from 
one of the finite- TV singularities of uj^ n > which survives as N — > 00 because it becomes asymptotically fiat, i.e., the 
determinant of the Hessian of W n ' goes to zero. This is the condition for a finite- TV singularity to survive given in 
However, preliminary results for n — 2 and d = 2, 3 [3(j do not support this scenario. 
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